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The Majorana fermion, a particle which act as its own an¬ 
tiparticles, was suggested by Majorana in 1937. A historic 
review of this problem can be found in lit]. The unambigu¬ 
ous observation of Majorana fermions remains an outstanding 
goal. No fundamental particles are known to be Majorana 
fermions, although there are speculations that the neutrino 
may be one. There are also theoretical suggestions that Majo¬ 
rana fermions may comprise a large fraction of cosmic Dark 
Matter fl. It has also been predicted that certain topological 
media, which are gapped in the bulk, may have topologically 
protected surface excitations with the properties of Majorana 
fermions fUj. This class of topological media includes topo¬ 
logical insulators such as superfluid ^He-B, thin films of ^He- 
A, vortex lines in superfluids and superconductors 0], thin su¬ 
perconducting wires and the quantum vacuum of the Standard 
Model of the Universe 0. Graphene may also belong to this 
class if a broad enough energy gap may appear due to spin- 
orbit interaction (or due to spontaneously broken symmetry) 
0. In this letter we report the first direct observation of gap¬ 
less Majorana quasiparticles which appear as Andreev bound 
states on the surface of superfluid ^He-B 0 . We have made 
precise measurements of the heat capacity of superfluid ^He-B 
at the extreme low-temperatures limit. We are able to separate 
the heat capacity contributions of bulk Bogolyubov quasipar¬ 
ticles and the surface Majorana quasipartilces by their differ¬ 
ent temperature dependencies. We have found that at 0.11 mK 
the Majorana fermions contribute about half of the bolometer 
heat capacity under the conditions of our experiments. 


There are several scenarios for the experimental search for 
Majorana signatures iffl. Many investigations have fo¬ 
cused on thin supeconducting wires. Here, the region for Ma¬ 
jorana formation is situated near the ends of wires, meaning 
that the Majoranas would have exactly zero energy and would 
be unable to move. The question is: how would it be possi¬ 
ble to observe a not-moving zero-energy state, a difficult task. 
There have been attempts to confirm the existence of such Ma¬ 
joranas indirectly by the observation of interference of Majo¬ 
rana states on each side of the wire through a supercondacting 
bridge 0 . However, a better choice is offered by the topolog¬ 
ical insulator, the time-reversal invariant B-phase of superfluid 
^He. The spin-triplet superfluid supports the existence of Ma¬ 
jorana quasiparticles 1 12 , El. This follows from the particle- 
hole symmetry of the Bogoliubov quasiparticles, = J-e- 
Consequently at zero energy they satisfy the Majorana condi¬ 
tion Yq = Yq- Zero-energy quasiparticle bound states appear 


in ^He-B when the underlying quasiparticle potential changes 
sign. This conditions is satisfied at edges, at surfaces and on 
vortices with odd integer winding numbers, where the quasi¬ 
particles form zero-energy Andreev bound states (a unique 
quasiparticle state in which Andreev reflection li, 14] plays 
a fundamental role). 

The superfluid energy gap of ^He-B is suppressed near 
walls over a distance corresponding to the superfluid coher¬ 
ence length which is about 80 nm at zero bar. In contrast to 
superconducting wires, the Majorana quasiparticles in ^He-B 
are able to move along the surface of the sample and con¬ 
sequently have a kinetic energy. They have a linear disper¬ 
sion relation with a zero gap, referred to as the Majorana cone 
ilEi. Because of this specific dispersion relation, the Majo¬ 
rana quasiparticles have a non-zero energy and consequently 
have a heat capacity. A proposed experimental method for 
observing the additional Majorana heat capacity can be found 
in HI, and the first attempt of its realization in IT^ . In the 
latter case the experiments were performed in superfluid ^He 
in a micron powder heat exchanger at a temperatures near Tc. 
Unfortunately, it is impossible to apply the Majorana theory 
quantitatively to the results of these experiments. 

The Majorana heat capacity obeys a power law with tem¬ 
perature, with a quadratic dependence on T for the specular 
scattering of quasiparticles. The Bogolyubov quasiparticles 
in the bulk ^He have a gap and consequently the bulk heat ca¬ 
pacity falls exponentially with reducing temperature. At some 
low temperature the heat capacity of the Bogolyubov quasi¬ 
particles becomes lower than that of the Majoranas one. This 
temperature depends on the ratio between the volume and the 
surface of the ^He sample in the bolometer. In the experiments 
described in this letter, we have observed a deviation from an 
exponential law for the quasiparticles heat capacity in bolome¬ 
ter, which corresponds well to an additional heat capacity due 
to zero-gapped Majorana quasiparticles. The very good agree¬ 
ment with theoretical predictions supports our claim of a di¬ 
rect observation of zero-gapped Majorana quasiparticles. 

Fig.l shows schematically the conditions of our experi¬ 
ments. The experiments are made at zero pressure and at tem¬ 
peratures down to 0.11 mK. The superfluid ^He-B energy gap 
at zero pressure is about 2 mK. Thus the quasiparticles gas is 
very dilute and consequently its heat capacity is very small. 
Near the wall the energy gap falls to zero. This region allows 
the existence of Majorana quasiparticles with zero energy gap. 
These quasiparticles can move in the 2D space along the sur- 
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FIG. 1: The experimental scenario. On the x axis we plot the dis¬ 
tance from the ^He boundary in coherence length units. The energy 
gap is shown by the solid line and temperature by the dashed line. 
Quasiparticles with energies above the gap can move in the bulk 
^He. States with zero gap, the Majorana particles, only exist near 
the wall, but can move transversely along the walls with a kinetic 
energy which contributes an additional heat capacity. 


face. The heat capacity of bulk Bogolyubov quasiparticles and 
Majorana quasiparticles can be estimated as followings: 

The Bogolyubov quasiparticles heat capacity in superfluid 
^He-B falls with temperature as OlTlIlSll : 


Cbuik (A)3/2exp (_ A) 


( 1 ) 


where Pp is the fermi momentum, A the superfluid gap ~ 2kTc 
and V is the volume of the sample. 

Owing to the ze ro g ap the Majorana heat capacity follows 
the power law lElIil: 
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where ^ is the ^He-B coherence length and A the surface area 
of the sample. 

The ratio of these heat capacities, including the numerical 
factors, reads: 
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where A is a geometrical parameter characteristic of the ex¬ 
perimental cell, the ratio between the volume and the surface 
area. The function F, introduced by this equation, is very use¬ 
ful since it indicates the temperature where the Majorana heat 
capacity begins to become significant. In a current experi¬ 
ments the parameter Ais about 10^^. The crossover when 
the heat capacities of the Majorana and Bogolyubov quasi¬ 
particles are equal should occur when F = 10^, that is about 
0.105 mK 12^. which corresponds well with our experimental 
results. 

We have used a bolometer which consists of a closed cop¬ 
per box with a small orifice lE^ . The box has the form of a 


cylinder of 6 mm diameter and 5 mm height. The diameter 
of the orifice is about 0.2 mm. The top and bottom plates of 
the cell are made from copper foils and the cylinder is turned 
copper. The volume is about 0.13 cm^. The bolometer is situ¬ 
ated inside the chamber of a nuclear demagnetization refriger¬ 
ator filled with superfluid ^He at extremely low temperatures. 
The only thermal reservoir in the bulk superfluid ^He consists 
of the thermal gas of Bogolyubov quasiparticles. The tem¬ 
perature inside the bolometer is determined by the balance of 
heat leak in and cooling by the flow of quasiparticles flow¬ 
ing out of the orifice. Because of the poor Kapitza resistance, 
the thermal conduction between the ^He and the walls of the 
bolometer can be taken as zero. 

To calculate the heat capacity we measure the thermal re¬ 
sponse of the liquid in the bolometer to a calibrated heat pulse. 
To measure the temperature (density of Bogolyubov quasipar¬ 
ticles) we have used the vibrating wire resonator (VWR) tech¬ 
niques first described in 112111 . The broadening of the VWR 
resonance is determined by the dissipative interaction with the 
Bogolyubov quasiparticles. (For a description of the VWR 
see Methods.) There are two different temperature scales for 
ther frequency width of a VWR as a function of temperature, 
the Lancaster scale 12 If and the Grenoble scale 12211 . Here 
we have used the Lancaster scale which fits our experimental 
data better. After an energy deposition the quasiparticle den¬ 
sity rises rapidly and then decreases with a time scale, of a 
few seconds, as shown in Fig. 2. The time constant Xb is deter¬ 
mined by the ratio of the volume of the bolometer and the ar- 
era of the orifice. If we take into account only the Bogolyubov 
quasiparticles, then Xb changes only slowly with cooling since 
the heat capacity and energy flow through the orifice are both 
proportional to the quasiparticle density. However, the time 
constant changes as soon as the additional heat capacity aris¬ 
ing from the Majorana quasiparticles appears. By observing 
the temperature-dependence of Xb we have made the first con¬ 
firmation of the existence of Majorana quasiparticles in super¬ 
fluid ^ He S El. 

The quasiparticle density and temperature increase almost 
instantaneously after a heating pulse. The measured VWR 
broadening increases with some time delay due to the quality 
factor of the VWR. The response time Xr is inversely propor¬ 
tional to the baseline width of the VWR (Wq); Xr = 1/{nWo). 
Consequently, the time dependence of measured VWR re¬ 
sponse on the heating pulse at to can be written 0: 


W{t)=Wo+A 
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where A is the fitting parameter, the broadening of VWR im¬ 
mediately after the heat pulse, for the conditions of zero time 
delay of VWR thermometer. We have fitted our data using 
this function as shown in Fig. 2 and tabulated the values of A. 
We note that the treatment of the data in this letter is differ¬ 


ent from that in Il24l] . Here we have calculated the value of A 
separately for each event. 

Direct heating of superfluid ^He by an immersed heater is 
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FIG. 2: The evolution of the VWR width after a heating pulse. The 
red line is a fit of the experimental points by Eq. (8). The blue line is 
the corresponding quasiparticle density obtained by Eq. (7). H is the 
actual measured peak width increase and A the corresponding value 
assuming instantaneous VWR response. 
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FIG. 3: The amplitude of heater VWR oscillation just after the RF 
pulse (stars). The straight line shows the linear regime of quasiparti¬ 
cle excitations, as discussed in the text. Bullets shows the additional 
broadening of the thermometry VWR A after heating pulses of dif¬ 
ferent energies. The baseline width is about 0.39 Hz. 


impossible owing to the very large Kapitza resistance between 
heater and liquid. Heat can be injected directly into the super¬ 
fluid by exciting a VWR beyond the critical velocity for pair 
breaking as proposed in 11251] . One should keep in mind that 
the electrically injected power is firstly transformed into me¬ 
chanical energy of the VWR before being transferred to the 
liquid. By pair breaking the VWR motion creates the quasi¬ 
particles directly in the superfluid ^He each with a momen¬ 
tum of ^ Pp. By integrating the voltage and current through 
the VWR we can calculate the deposited energy. The question 
is: how much of this energy is transformed to quasiparticles? 
In our previous Dark Matter detector investigations 11241] we 


FIG. 4: The calibration factor for a range of temperatures. The 
dashed line represents the variation expected from the Bogolyubov 
quasiparticles alone. The solid line takes into account the additional 
heat capacity of the Majorana quasiparticles. 


T, pK 


100 no 120 130 140 150 160 170 180 



FIG. 5: The heat capacity ratio C/Chuik calculated from data shown 
in Fig.4. The line equal to 1 correspond to Bogoluybov QPs. The 
curved dashed line is that for the Majorana quasiparticles alone, and 
the solid line corresponds to a heat capacity of a both components. 


surprisingly found a decrease of the sensitivity on cooling be¬ 
low 200 jj.K. To explain this feature we suggested that part of 
VWR heater energy is lost through the intrinsic damping of 
the VWR. We suggested that a fraction of the energy W,„f/W q 
was lost due to VWR heating. Indeed, the parameter Wim was 
a fitting parameter. In our current experiments we have found 
that this mechanism cannot explain the increasing of time con¬ 
stant Tp, we observed earlier 1201] . Furthermore, we have ana¬ 
lyzed the oscillation amplitude of the heating VWR just after 
the heating pulse. We find a clear saturation of the ampli¬ 
tude. In other words, the VWR reaches the critical velocity for 
cooper pair breaking as shown in Fig. 3. Thus the mechanical 
energy transfer to the quasiparticles is very effective and we 
should assign all the dissipated VWR energy as contributing 
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to the calibration pulses. The decrease of sensitivity, observed 
earlier, cannot be explained by mechanical loses, but rather by 
the additional heat capacity arising from Majorana particles. 
In Fig. 3 the calibration of the additional VWR broadening af¬ 
ter the heating VWR excitation over a range of energies is also 
shown. By these measurements at different base temperature 
we are able to calculate the calibration factor of the bolome¬ 
ter, as shown in Fig. 4. The dashed line shows the calibra¬ 
tion factor as a function of temperature (VWR damping) for 
the Bogolyubov quasiparticle heat capacity alone. However, 
we see a significant deviation of the experimental data from 
this line, but if we also add the heat capacity of the Majorana 
quasiparticle (Eq. 2) we hnd good agreement between the ex¬ 
perimental data and the theoretical curve. In Fig. 5 the ratio of 
the experimental data and the Bogolyubov quasiparticle heat 
capacity is shown. We know the volume of our bolometer, 
but the surface of the walls has to be estimated taking into ac¬ 
count the surface roughness. The walls of the bolometer are 
of turned copper (cylindrical body) and copper foils (ends). In 
both cases the roughness of the walls may contribute an effec¬ 
tive area 5-10 times larger than the geometrical surface area. 
Electron microscope analysis of the surface also shows high 
roughness. We have used the surface area of the walls as a 
fitting parameter and found that it to be about 10 time bigger 
than the geometrical surface, in a reasonable agreement with 
the apparent roughness. 


We should comment on the influence on our results of a 
possible layer of solid ^He on the surfaces of the bolometer. 
In our case the surface was covered by ^He during the con¬ 
densation of the liquid ^He into the cell. Owing to the small 
surface-to-volume ratio it is enough to use the 0.1% of ^He 
to completely prevent any solid ^He layer to form inside the 
bolometer. It was shown in 11271] that in the presence of frac¬ 
tion of solid ^He one should observed the two bolometer time 
constants arising from the the imperfect thermal contact be¬ 
tween the solid and liquid ^He. We have not seen any such 
signature of solid ^He in our current experiments. 


We thus can draw the conclusion that the existence of Ma¬ 
jorana fermions is confirmed in this system by the very direct 
method of measuring the added contribution of the Majorana 
heat capacity. 

METHODS 


The VWR comprises a superconducting NbTi wire hlament 
bent into an approximately semi-circular shape of a few mm 
diameter, with both ends hrmly fixed lE^ . The VWR is situ¬ 
ated in a small magnetic field and driven by the Lorentz force 
generated by an AC current passed around the loop at close 
to the mechanical resonance frequency, oscillating perpendic¬ 
ularly to its main plane with an rms velocity v. The motion 
is damped by frictional forces of total amplitude F'(v) aris¬ 
ing mainly from momentum transfer to the quasiparticles of 
the surrounding superfluid with a magnitude proportional to 
the quasiparticle density. The quasiparticles density depends 
exponentially on the temperature so that the temperature de¬ 


pendence of the VWR damping can be written as Kill : 

A 


W {T) = a exp ( — 


kT 


( 5 ) 


where a is constant determined by the microscopic properties 
of the liquid. The rapid exponential function ensures that the 
VWR provides a very sensitive thermometer. Since it has been 
found that the damping depends slightly on the wire velocity 


1241281], we have used very small amplitudes of VWR exci¬ 
tation, which we have adjusted with temperature to ensure a 
constant amplitude of the VWR resonator. 
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